qr_rank1_update

public interface qr_rank1_update

Module Procedures

private pure subroutine qr_rank1_update_dbl(q, r, u, v)

Computes the rank-1 update to an M-by-N QR factored matrix $A$ where $M \ge N$ , $A = Q R$ , and $A_1 = A + \vec{u} \vec{v}^T$ such that $A_1 = Q_1 R_1$ .

Arguments

Type	Intent	Attributes		Name
real(kind=real64),	intent(inout),	dimension(:,:)	::	q	On input, the original M-by-K orthogonal matrix $Q$ . On output, the updated matrix $Q_1$ .
real(kind=real64),	intent(inout),	dimension(:,:)	::	r	On input, the M-by-N matrix $R$ . On output, the updated matrix $R_1$ .
real(kind=real64),	intent(in),	dimension(:)	::	u	The M-element $\vec{u}$ update vector.
real(kind=real64),	intent(in),	dimension(:)	::	v	The N-element $\vec{v}$ update vector.

private pure subroutine qr_rank1_update_cmplx(q, r, u, v)

Computes the rank-1 update to an M-by-N QR factored matrix $A$ where $M \ge N$ , $A = Q R$ , and $A_1 = A + \vec{u} \vec{v}^H$ such that $A_1 = Q_1 R_1$ .

Arguments

Type	Intent	Attributes		Name
complex(kind=real64),	intent(inout),	dimension(:,:)	::	q	On input, the original M-by-K orthogonal matrix $Q$ . On output, the updated matrix $Q_1$ .
complex(kind=real64),	intent(inout),	dimension(:,:)	::	r	On input, the M-by-N matrix $R$ . On output, the updated matrix $R_1$ .
complex(kind=real64),	intent(in),	dimension(:)	::	u	On input, the M-element $\vec{u}$ update vector. On output, the original content of the array is overwritten.
complex(kind=real64),	intent(in),	dimension(:)	::	v	On input, the N-element $\vec{v}$ update vector. On output, the original content of the array is overwritten.

qr_rank1_update Interface

Contents

Module Procedures

public interface qr_rank1_update

Module Procedures

private pure subroutine qr_rank1_update_dbl(q, r, u, v)

Arguments

private pure subroutine qr_rank1_update_cmplx(q, r, u, v)

Arguments